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Abstract 

New data on proton and pion production in p+C interactions from the CERN PS and SPS 
accelerators are used in conjunction with other available data sets to perform a comprehen- 
sive survey of backward hadronic cross sections. This survey covers the complete backward 
hemisphere in the range of lab angles from 10 to 180 degrees, from 0.2 to 1.4 GeV/c in 
lab momentum and from 1 to 400 GeV/c in projectile momentum. Using the constraints 
of continuity and smoothness of the angular, momentum and energy dependences a con- 
sistent description of the inclusive cross sections is established which allows the control of 
the internal consistency of the nineteen available data sets. 
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1 Introduction 



An impressive amount of data on backward hadron production in p+C interactions has 
been collected over the past four decades. A literature survey reveals no less than 19 experiments 
which have contributed a total amount of more than 3500 data points covering wide areas in 
projectile momentum, lab angle and lab momentum. 

Looking at the physics motivation and at the distribution in time of these efforts, two 
distinct classes of experimental approaches become evident. 15 experiments cluster in a first 
period during the two decades between 1970 and 1990. All these measurements have been 
motivated by the nuclear part of proton-nucleus collisions, in particular by the width of the 
momentum distributions in the nuclear rest system which reach far beyond the narrow limits 
expected from nuclear binding alone. These studies have ceased in the late 1980's with the 
advent of relativistic heavy ion collisions and their promise of "new" phenomena beyond the 
realm of classic nuclear physics. 

A second class of very recent measurements has appeared and is being pursued after 
the turn of the century, with publications starting about 2008. Here the motivation is totally 
different. It is driven by the necessity of obtaining hadronic reference data for the study of 
systematic effects in cosmic ray and neutrino physics, in particular concerning atmospheric and 
long base line experiments as well as eventual novel neutrino factories. The main aim of these 
studies is the comparison to and the improvement of hadronic production models - models 
which are to be considered as multi-parameter descriptions of the non-calculable sector of the 
strong interaction, with very limited predictive power. 

This new and exclusive aim has led to the strange situation that if all recent publications 
contain detailed comparisons to available production models, no comparison to existing data is 
attempted. It remains therefore unclear how these new results compare to the wealth of already 
available data and whether they in fact may over-ride and replace the existing results. 

In this environment the studies conducted since 15 years by the NA49 experiment at the 
CERN SPS have a completely different aim. Here it is attempted to trace a model-independent 
way from the basic hadron-nucleon interaction via hadron-nucleus to nucleus-nucleus colli- 
sions. This aim needs precision data from a large variety of projectile and target combinations 
as well as a maximum phase space coverage. As the acceptance of the NA49 detector is limited 
to lab angles below 45 degrees, it is indicated to use existing backward data in the SPS energy 
range in order to extend the acceptance coverage for the asymmetric proton-nucleus interac- 
tions. This requires a careful study of the dependence on cms energy and of the reliability of the 
results to be used. 

In the course of this work it appeared useful and even mandatory to provide a survey 
of all available data over the full scale of interaction energies, the more so as no overview of 
the experimental situation is available to date. This means that the present study deals with 
projectile momenta from 1 to 400 GeV/c, for a lab angle range from 10 to 180 degrees, and for 
lab momenta from 0.2 to 1.2 GeV/c. 



2 Variables and kinematics 

Most available data have been obtained as a function of the lab momentum piab (or kinetic 
energy Tiab) at constant lab angle ©lab. In this publication all given yields are transformed to the 
double differential invariant cross section 

/(Plab, Qlab) = . ^ • 
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In this context the term "backward" needs a precise definition. One possibility would be 
to define as "backward" the region of lab angles ©lab > 90 degrees. The present paper uses 
instead a definition which refers to the cms frame with the basic variables Feynman xp and 
transverse momentum p-r, defining as "backward" the particle yields dXxp < 0. This allows a 
clear separation of the projectile fragmentation region at positive xp with a limited feed-over 
into negative xp and the target fragmentation region at negative xp with a limited feed-over 
into positive xp. M the same time the notion of "kinematic limit" in participant fragmentation 
is clearly brought out dX xp = ±. 1 and the contributions from intranuclear cascading may be 
clearly visualized and eventually separated. 

The correlation between the two pairs of variables is presented in Fig. [T] which shows 
lines of constant piab and ©lab in the coordinate frame of a;^ and pt for protons and pions for the 
two values of projectile momentum at 158 and 3 GeV/c which are representative of the typical 
range of interaction energies discussed in this paper. 

Several comments are due in this context. The definition of Feynman x p has been modi- 
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Figure 1: Lines of constant piab and 0iab in the cms frame spanned by Feynman xp and pp for 
protons and pions at two different projectile momenta, a) protons at 158 GeV/c, b) protons at 
3 GeV/c, c) pions at 158 GeV/c and d) pions at 3 GeV/c 
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fled from the standard one, 
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Xp = 




with rup the proton mass. This takes care of baryon number conservation and regularises the 
kinematic borders at low interaction energies. The s dependence in Fig.[T]is small to negligible 
for lab angles above about 50 degrees both for pions and protons but becomes noticeable at 
small ©lab. If at SPS energy the full range of lab momenta up to 1.4 GeV/c and angles above 
10 degrees is confined to the backward region both for protons and pions, the coverage for pions 
extends to positive xp at low lab angles and low beam momenta. 

Another remark concerns the overlap between target fragmentation and nuclear cascad- 
ing. For protons, at all lab angles above about 70 degrees the kinematic limit for fragmentation 
of a target nucleon at rest in the lab system is exceeded. For pions on the other hand this is not 
the case as their xp value for p^^b = is at 



This means that over the full range of lab angles and up to large piab values the contribution 
from target participants mixes with the nuclear component. The separation of the two processes 
therefore becomes an important task, see Sect. [10] of this paper. 

A last remark is due to the limits of experimental coverage. All existing experiments run 
out of statistics at cross section levels of about 10 /ib, that is about 4 orders of magnitude below 
the maximum yields. As visible from the momentum ranges indicated in Tables 1 and 2, this 
corresponds to a typical upper momentum cut-off in the region of 1 GeV/c. 

3 The Experimental Situation 

The backward phase space coverage in p-i-C interactions is surprisingly complete if com- 
pared with the forward direction and even with the available data in the elementary p+p colli- 
sions. This is apparent from the list of experiments given in Tables [T] and |2] with their ranges in 
beam momentum, lab angle, and lab momentum. Although some effort has been spent to pick 
up all published results, this list is not claimed to be exhaustive as some results given as "private 
communication", in conference proceedings or unpublished internal reports might have escaped 
attention. 

For secondary protons. Table [B the important amount of low energy n-i-C data by Franz 
et al. [10 1 has been added to the survey as the isospin factors for the transformation into p+C 
results have been studied and determined with some precision, see Sect. [51 

For secondary pions. Table [2l the situation is somewhat complicated by the fact that two 
independent sets of results have been published by the HARP-CDP [4] and the HARP |[T5l 
groups, based on identical input data obtained with the same detector. An attempt to clarify this 
partially contradictory situation is presented in Sect. |9.3| of this paper. 

Unfortunately, no commonly agreed scale in the three basic variables ©lab, Puh and Pbeam 
of the double-differential cross sections has been defined by the different collaborations pro- 
viding the data contained in Tables [U and [21 This leads to the fact that not a single couple out 
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interaction 


Experiment 


projectile momentum 


lab angle coverage 


Plab coverage 


number of 


errors [%] 






(GeV/c) 


(degrees) 


(GeV/c) 


data points 


(cstat) 


(cTsyst) 




Bayukov [T] 


400 


70, 90, 118, 137, 160 


0.4-1.3 


35 


6 


20 




NA49 


158 


10, 20, 30, 40 


0.3-1.6 


40 


7 


5 


p+C 


Belyaev |3| 


17, 23,28,34 ,41,49, 56 


159 


0.3-1.2 


125 


5 


15 


HARP-CDP tU 


3,5, 8, 12, 15 


25,35,45, 55,67, 82, 97, 112 


0.45-1.5 


202 


4 


6 




Burgov (5) 


2.2, 6.0, 8.5 


162 


0.35-0.85 


36 


15 


5 




Bayukov fSl 


1.87,4.5,6.57 


137 


0.3-1.1 


55 


10 


20 




Geaga [7 J 


1.8,2.9,5.8 


180 


0.3-1.0 


50 


17 


15 




Frankel (8) 


1.22 


180 


0.45-0.8 


6 


7 






Komarov f9\ 


1.27 


105. 115, 122, 130, 140, 150, 160 


0.34-0.54 


-200 


8 


15 


n+C 


Franz 1 10 | 


0.84, 0.99, 1.15 


51, 61, 73, 81, 98, 120, 140, 149, 160 


0.3-0.8 


553 


5 


10 



Table 1 : Data sets for proton production in p+C and n+C collisions from seven experiments 
giving the ranges covered in projectile momentum, lab angle, and lab momentum, the number 
of measured data points and errors 



Experiment 


projectile momentum 


lab angle coverage 


piab coverage 


number of 


errors [%] 




(GeV/c) 


(degrees) 


(GeV/c) 


data points 


{cstat 




Nikiforov | 11 1 


400 


70, 90, 118, 137, 160 


0.2-1.3 


59 


12 




NA49 ri2j 


158 


5, 10, 15. 20, 25, 30, 35, 40, 45 


0.1-1.2 


174 


5 


4 


Belyaev [Uj 


17,22, 28, 34,41,47,57 


159 


0.25-1.0 


218 


4 


15 


Abgrall (U) 


31 


0.6-22.3 


0.2-18 


624 


6 


7 


HARP-CDP (4) 


3, 5, 8, 12, 15 


25, 35,45,55,67, 82, 97, 112 


0.2-1.6 


829 


6 


8 


HARP [15] 


3,5, 8, 12 


25, 37, 48, 61, 72, 83, 95, 106,117 


0.125-0.75 


605 




12 


Burgov [l6] 


2.2, 6.0, 8.5 


162 


0.25-0.6 


29 




20 


Baldin Qt] 


6.0, 8.4 


180 


0.2-1.25 


45 




10 


Cochran fT8l 


1.38 


15, 20, 30, 45, 60, 75, 90, 105, 120, 135, 150 


0.1-0.7 


199 


3 


12 


Crawford |19 | 


1.20 


22.5, 45, 60, 90, 135 


0.1-0.4 


50 


8 


7 



Table 2: Data sets for pion production in p+C collisions from seven experiments giving the 
ranges covered in projectile momentum, lab angle, and lab momentum, the number of measured 
data points and errors 



of the more than 3500 data points contained in these Tables may be directly compared. The 
application of an interpolation scheme as described in Sect.|4]is therefore an absolute necessity. 
Ideally the thus obtained interpolated cross sections would form an internally consistent sam- 
ple of results which would be coherent within the given experimental errors. As will become 
apparent in the following data comparison, this assumption is surprisingly well fulfilled for the 
majority of the experiments. Only four of the 20 quoted groups of results fall significantly out of 
this comparison; those will be discussed in Sect.|9]of this paper. In this sense the overall survey 
of the backward proton and pion production results in a powerful constraint for the comparison 
with any new data sample. 

4 Data comparison 

As stated above the main problem in bringing the wealth of available data into a consistent 
picture is given by the generally disparate position in phase space and interaction energy of the 
different experiments. The triplet of lab variables given by the beam momentum Pbeam. the lab 
momentum piab and the lab angle ©lab has been used for the establishment of the following 
interpolation scheme. In addition and of course, the statistical and systematic errors have to be 
taken into account in the data comparison. 
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4.1 Errors 



The last columns of Tables \T\ and [21 contain some information about the statistical and 
systematic errors of the different experiments. The given numbers are to be regarded as mean 
values excluding some upward tails as they are inevitable at the limits of the covered phase 
space in particular for the statistical uncertainties. In some cases only rudimentary information 
about the systematic errors is available or the systematic and statistical errors are even combined 
into one quantity. In the latter cases these values are given in between the respective columns of 
Tabled 

Inspection of these approximate error levels reveals a rather broad band of uncertainties 
ranging from about 4% to about 20%, the latter limit being generally defined by overall normal- 
ization errors. The presence of extensive data sets well below the 10% range of both statistical 
and systematic errors gives however some hope that a resulting overall consistency on this level 
might become attainable by the extensive use of data interpolation. 

The term "interpolation" is to be regarded in this context as a smooth interconnection of 
the data points in any of the three phase space variables defined above. This interconnection is 
generally done by eyeball fits which offer, within the error limits shown above, sufficient accu- 
racy. If the distributions in 6iab and interaction energy are anyway not describable by straight- 
forward arithmetic parametrization, thepiab dependences are, as discussed in Sect. l4.4l below. in 
a majority of cases approximately exponential. In these cases exponential fits have been used if 
applicable. 

As additional constraint physics asks of course for smoothness and continuity in all three 
variables simultaneously. Therefore the resulting overall data interpolation has to attempt a 
three-dimensional consistency. 

If the data interpolation helps, by the inter-correlation of data points, to reduce the local 
statistical fluctuations, it does of course not reduce the systematic uncertainties. It is rather 
on the level of systematic deviations that the consistency of different experimental results is 
to be judged. It will become apparent from the detailed discussion described below that the 
majority of the quoted experiments allows for the establishment of a surprisingly consistent 
overall description in all three variables. 

4.2 Dependence on cms energy s 

As the data discussed here span an extremely wide range of cms energy from close to 
production threshold to the upper range of Fermilab energies, a suitable compression of the en- 
ergy scale has been introduced in order to be able to present the results in a close-to-equidistant 
fashion against energy. The form chosen here is the variable 1 / y/s. This choice is suggested by 
the considerable amount of work invested in studying the approach of hadronic cross sections to 
the scaling limit at high energy in the 1970's [20J. In fact the Regge parametrization suggested 
a smooth dependence of the cross sections as s~", with a = 0.25-0.5 depending on the choice 
of trajectories involved. Such behaviour was indeed found experimentally. In the present study 
the cross sections turn out to have only a mild 1 / y/s dependence for -^i > 5 GeV, a dependence 
which is however different for pions and protons. This dependence is strongly modified below 
^/s ~ 2.5 GeV due to threshold effects. 

4.3 Angular dependence 

A convenient and often used scale for the lab angle dependence is given by cos (9iab). 
This scale has the advantage of producing shapes that are again to zero order exponential. Of 
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course, continuity through 6iab = 180 degrees imposes an approach to 180 degrees with tangent 
zero. As the data samples are generally not measured at common values of 9iab, a fixed grid of 
angles has been defined based on the 6iab values of the HARP-CDP experiment [4] dominating 
the range from 25 to 112 degrees. Measured values down to 10 degrees and in the higher angular 
range at 137, 160, and 180 degrees have been added. Measurements not corresponding to these 
grid values are interpolated using the cos (6iab) distributions specified below. 

4.4 Lab momentum dependence 

All data discussed here have been transformed into invariant cross sections ©• This fa- 
cilitates the presentation in different coordinate systems and eliminates the trivial approach of 
the phase space element to zero with decreasing momentum. In addition, most of the invariant 
p\ah distributions are close to exponential within the measured pi^h range. There are notable de- 
viations mostly at low momentum and in the lower (higher) range of lab angles for pions and 
protons, respectively, as well as in the approach to threshold. In these cases an eyeball fit has 
been used which can be reliably performed within the error margins indicated above. 

At low lab momenta physics requires a deviation from the exponential shape as the in- 
variant cross sections must approach piab = with tangent zero. This limit appears in general 
at piab < 0.2 GeV/c for pions and piab < 0.5 GeV/c for protons. The data presented here fall 
practically all above these momentum limits. Only the HARP experiment [fT5ll gives results at 
Piab = 0.125 GeV/c for pions where indeed a substantial deviation from the exponential shape 
is visible. This is shown in Fig. [2] where the deviation from exponential fits at this piab is given 
in percent for all angles and beam momenta together with similar deviations observed in p+p 
interactions \2U . 

p^ = 0.125 GeV/c; 
NA49 [21] J 



\ 

\ 

20 40 60 0.1 0.2 0.3 

A [%] Xp 

Figure 2: a) Distribution of the deviation of the data points at piab = 0.125 GeV/c from the 
exponential fits for 7r+ and tt^ at all angles and beam momenta, b) Deviation of 7r+ cross 
sections at pr = 0.125 GeV/c from exponential fits to the higher pr region in p-i-p interactions 
as a function of xp 

A number of examples of momentum distributions for protons and pions is given in the 
following Figs. [3] and [5] which show the invariant cross sections as a function of piab and the 
corresponding exponential fits 
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/(Plab, 6lab, Pbeam) = ^(©lab, Pbeam) * exp(-piab/ B{Qi^b, Pbeam)) (5) 

which are, whenever necessary, supplemented by hand interpolations into the non-exponential 
regions. 

A first group of distributions in the medium angular range at 45 and 97 degrees is pre- 
sented in Fig. [3] for the HARP-CDP data concerning protons and pions, including exponential 
fits. 
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Figure 3: Invariant cross sections for protons, vr"*" and n~ as a function of piab at 0iab = 45 and 
97 degrees. Full lines: exponential fits. Broken line: hand-interpolation into the non-exponential 
region 



Evidently the exponential shape is within errors in general a good approximation to the 
momentum dependence. More quantitative information is contained in the normalized residual 
distributions of the data points, 

rnorm = A/cr (6) 

where A is the difference between data and fit and sigma the statistical error of the given data 
point. Should the fit describe the physics and should systematic effects be negligible, the dis- 
tribution of rnorm IS cxpcctcd to bc Gaussian with rms equal to unity. The rnorm distributions are 
given for the totality of the HARP-CDP data in Fig. gl 

These distributions are well described by centred Gaussians. The resulting rms values are 
however somewhat bigger than one signalling systematic experimental effects or a deviation of 
physics from the simple exponential parametrization. In view of the statistical errors of 4% to 
6% given by HARP-CDP (Tables [T] and [2l) these deviations are on the level of a few percent 
which is below the error margin to be anyway expected from the present general data survey. 
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Figure 4: Normalized residual distributions for protons and pions for the complete set of beam 
momenta and angles of the HARP-CDP data with the exception of a few points at low angles 
and momenta which clearly exhibit non-exponential behaviour 



Further examples of piab distributions from other experiments are given in Fig. [5] for a 
selection of particle type, beam momenta and angles. 




Figure 5: Several examples of invariant cross sections as a function of piab for a variety of particle 
type, lab angle and beam momenta including exponential fits (full lines) and, when necessary, 
eyeball fits into the non-exponential regions of pi^b (broken lines). 

Again the basically exponential shape of these distributions is evident. Characterizing the 
exponential fits by their inverse slopes -B(0iab7Pbeam) a smooth and distinct dependence on lab 
angle and beam momentum becomes visible as shown in Fig.[6j 

Compared to the strong dependence of B on ©i^b which ranges from 0.3 to 0.05 GeV/c, 
the only modest dependence on Pbeam of ~0.03 GeV/c for beam momenta from 3 to 158 GeV/c 
is noticeable. 

Following the above data parametrization a generalized grid of piab values between 0.2 
and 1.2 GeV/c, in steps of 0.1 GeV/c, may now be established. Concerning the lower and upper 
limits of this grid, an extrapolation beyond the limits given by the experimental values has been 
performed in some cases. This extrapolation does not exceed the bin width of the respective 
data lists and is therefore defendable in view of the generally smooth, gentle and well-defined 
piab dependences. 
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Figure 6: Inverse slopes i?(6iab,Pbeam) as a function of 6iab for four beam momenta from 3 to 
158 GeV/c, a) for protons, b) for 7r+ and c) for vr^. The full lines are drawn to guide the eye 



4.5 Physics constraints 

In the absence of theoretical predictability in the soft sector of the strong interaction, any 
attempt at bringing a multitude of experimental results into a common and consistent picture 
has to rely on a minimal set of model-independent physics constraints. In fact a "democratic" 
averaging of eventually contradictory data sets would only add confusion instead of clarity. 

4.5.1 Continuity 

Two examples of the continuity constraint have already been mentioned above: invariant 
Piah distributions have to approach zero momentum horizontally that is with tangent zero. The 
same is true for angular distributions in their approach to 180 degrees. 

4.5.2 Smoothness 

It is a matter of experimental experience in the realm of soft hadronic interactions that 
in general distributions in any kind of kinematic variable tend to be "smooth" in the sense 
of absence of abrupt local upwards or downwards variations. The widespread use of simple 
algebraic parametrizations has its origin in this fact, specifically in the absence of local maxima 
and minima, with the eventual exception of threshold behaviour of which some examples will 
become visible below. 



4.5.3 Charge conservation and isospin symmetry 

Charge conservation has of course to be fulfilled by any type of experimental result. This 
means for instance that for the interaction of a positively charged projectile (proton) with an 
isoscalar nucleus (Carbon) the /n^ ratio has to be greater or equal to unity over the full 
phase space invoking isospin symmetry (and of course the experience from a wide range of 
experimental results). The presence of data with /n^ < 1 therefore immediately indicates 
experimental problems. The inspection of tt^/tt^ ratios has the further advantage that a large 
part of the systematic uncertainties, notably the overall normalization errors, cancel in this ratio. 
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4.5.4 Isospin rotation of secondary baryons and projectile 

It has been shown that in proton induced nuclear collisions the yields of the secondary 
protons and neutrons are related by a constant factor of about 2 which is in turn related to 
the ratio of the basic nucleon-nucleon interaction [|22ll . Similarly, when rotating the projectile 
isospin from proton to neutron, it has been predicted that the yield ratio of secondary protons 
from proton and neutron projectiles 

pp/n ^ f{P + C^ P') 

f{n + C^p) ^'^ 

should be constant and equal to 2.5 for light nuclei |[23ll . The extensive and precise low-energy 
data set of Franz et al. ifTOll from n+C interactions has therefore been included in the present 
survey. These data present a welcome extension of the 1 / y/s scale into the region 0.47 to 0.49 
which is not covered for most of the angular range with proton projectiles. As shown below, 
these data fit indeed very well, after re-normalization, into the general 1 / y/s dependence of 
secondary protons where the low energy data by Frankel et al. [8] and Komarov et al. \9j at 
angles between 112 and 180 degrees provide an independent control of the normalization. 



4.5.5 Establishing a consistent set of data 

With these constraints in mind, and having established the parametrization and interpola- 
tion of the piah distributions as discussed above, one may now proceed to the attempt at sorting 
the 19 available experiments into a consistent global data set. It would of course be rather sur- 
prising if all experiments would fit into this global picture within their respective error limits. 
In fact it turns out that this procedure establishes a very strong constraint for possible devia- 
tions, as a large majority of results is creating a perfectly consistent picture both for protons 
and for pions. Only four of the 19 data sets cannot be brought into consistency with all other 
experiments without gravely affecting and contradicting the above constraints. These data are 
not included in the following global interpolation scheme. They will be discussed separately in 
Sect.mbelow. 



5 The proton data 

5.1 1 / a/s dependence 

The invariant proton cross sections are shown in Fig.|7]as a function of 1/y/s for a grid 
of ten lab angles between 25 and 180 degrees and constant lab momenta between 0.3 and 

1.2 GeV/c. The interpolated data points in each panel are identified by symbols correspond- 
ing to the different experiments. 

The solid lines are eyeball interpolations through the data points. A first remark concern- 
ing this Figure concerns the smoothness and continuity of the 1 / y/s dependences. The achieved 
overall consistency of all data is rather impressive even if single points are deviating in some 
areas of phase space. The salient features of the physics contained in these plots may be sum- 
marized as follows: 

- A strong yield suppression between l/y/s ~ 0.45 and the elastic limit at l/y/s = 0.53 is 
evident. 

- The n-i-C data [10| are well consistent with the p-i-C results in the overlap regions; they 
define a broad maximum of the cross sections at l/i/i ~ 0.46 at medium angles and low 

Plah- 
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Figure 7: Invariant cross sections for protons in p+C collisions as a function of 1/ y/s at fixed piab 
and ©lab- The interpolated data points are indicated by symbols corresponding to the respective 
experiments in each panel. The solid lines represent the interpolation of the data 

- There is a well-defined asymptotic behaviour of the cross sections for 1 / y/s below about 
0.2 or beam momenta above about 12 GeV/c. 

- For the lower 6iab region and/or low piab the asymptotic region is approached from above. 
The latter point is reminiscent of the behaviour of the proton yields in p+p interactions, 

as shown in Fig.[8l 

Another feature of Fig.|7]is the systematic droop of the cross sections from HARP-CDP at 
their highest beam momentum of 15 GeV/c or l/y/s = 0.18, demonstrating the discriminative 
power of the approach. This decrease is quantified in Fig. |9] where the ratio between the 
measured invariant cross sections and the data interpolation is shown as a function of pi^b for 
the complete angular range from 25 to 97 degrees. Here deviations of up to 50% are visible. 

The abruptness of this decrease would necessitate a rather violent variation of the cross 
sections with increasing energy including a minimum between PS and SPS energies. A final 
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Figure 8: Invariant proton cross sections as a function of 1/ for p+p interactions at the four 
lab angles a) 25, b) 35, c) 45 and d) 55 degrees for piab values from 0.4 to 1 GeV/c. The data are 
interpolated from Blobel \,24i and NA49 [25J. The lines are drawn to guide the eye 
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Figure 9: Ratio between the interpolated invariant proton cross sections from HARP-CDP 
and the global interpolation as a function of pi^b for the angular range 25 < ©lab < 97 degrees 



clarification of this situation is given by the proton data from Serpukhov [3] which, although 
suffering from a different and independent problem, at least exclude such variations in the region 
between 17 and 67 GeV/c beam momentum, see Sect. l9.2l below. 
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5.2 cos(6'iab) dependence 

In addition to the description of the energy dependence, the global interpolation has of 
course also to result in a smooth and continuous verification of the angular dependence present- 
ing the third dimension of the present study. This constraint has to be fulfilled at any value of 

In a first example the situation at I/a/s = 0.05 is shown in Fig. \W\ This value lies in 
between the Fermilab [IJ and NA49 [2] data in the region of negligible s -dependence. It there- 
fore allows for the direct comparison of the two experiments in their respective angular regions 
which have no overlap. 



160° 137° 118° 90° 70° 40° 30°20°10° 




COS i&J 

Figure 10: Invariant proton cross sections at l/y/s = 0.05 as a function of cos(6'iab) combining 
the Fermilab and NA49 data for pi^b between 0.2 and 1.4 GeV/c. The global interpolation is 
shown as full lines. The measured cross sections in the angular ranges from 70 to 160 degrees |[T1 
and from 10 to 40 degrees [2 J are given on the vertical broken lines 

Several observations are in place here: 

- The two experimental results connect perfectly through the gap between the NA49 
(0\ab < 40 degrees) and the Fermilab (^^lab > 70 degrees) data. 

- There is at most a few percent variation of the cross sections between the angles of 160 
and 180 degrees taking into account the constraint of continuity through 180 degrees 
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discussed in Sect. |4] above. This allows the combination of results in this angular region 
as it is applied in the determination of the 1 / y/s dependence, Fig. U\ 

- The angular distributions are smooth and close to exponential in shape. In particular, no 
instability in the region around 90 degrees is visible where an eventual diffractive peak 
from target fragmentation would appear, see also El. 

Further angular distributions at four 1 / ^/s values between 0. 1 and 0.4 GeV^^ are given in 
Fig. [TT] In fact such distributions at arbitrary values of 1/ ^/s may be obtained from the global 
interpolation as it is presented in numerical form at the NA49 web page [|26ll . 




Figure 1 1 : Invariant proton cross sections as a function of cos(6'iab) for four values of 1/ y/s: a) 
0.1, b) 0.2, c) 0.3, d) 0.4 GeV"^ and for piab values between 0.4 and 1.2 GeV/c. The standard 
grid of 10 angles, Fig.|71 is indicated by the vertical broken lines 



Evidently the angular distributions maintain their smooth and continuous shape, specifi- 
cally through 90 degrees, at all interaction energies. With the approach to low beam momenta 
however, a progressive rounding of the shape towards higher lab angles manifests itself. 
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The data for positive pions 



The global interpolation of the 7r+ data is presented in this section in close analogy to the 
preceding section for protons. 

6.1 1 / a/s dependence 

The invariant 7r+ cross sections are shown in Fig.[12]as a function of 1/ y/s for the standard 
grid of ten lab angles between 25 and 1 80 degrees and for constant lab momenta between 0.2 and 

1.2 GeV/c. The interpolated data points in each panel are identified by symbols corresponding 
to the different experiments. 

The solid lines represent the global interpolation by eyeball fits of both the energy and the 
angular dependences. Again the l/^/s dependence is in general smooth and continuous, with 
an impressive overall consistency of all data with only few exceptions discussed below. There 
are some general trends to be pointed out: 

- At the lowest lab momentum, the pion cross sections are practically s-independent, with 
variations of only 10-20% in the range from 1 to 400 GeV/c beam momentum. 

- This fact suggests 7r+ production at low momentum transfer in the nuclear cascade. 

- For all lab momenta, the approach to high energies is very flat for I/a/s < 0.2 or beam 
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Figure 12: Invariant cross sections for 7r+ in p+C collisions as a function of 1/ y/s at fixed piab 
and 6'iab. The interpolated data points are indicated by symbols corresponding to the respective 
experiments in each panel. The solid lines represent the global data interpolation 
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momenta above 12 GeV/c. 
- The high energy cross sections are approached for all angles and beam momenta from 
below. 

There are two areas of deviation from the global interpolation which are both connected 
to the HARP-CDP data [4J. At their lowest angle of 25 degrees, the cross sections are system- 
atically low by up to a factor of two below piab ~ 0.5 GeV/c and 1/ ^/s above 0.2. This is in 
contradiction to the available low energy data from other experiments also shown in Fig. [121 
The second area concerns, as for the protons, the data at 15 GeV/c beam momentum where a 
characteristic pattern of deviations is visible: At low angles and low piab, the data tend to over- 
shoot the interpolation, whereas at angles above 45 degrees a progressive droop with increasing 
lab momentum is evident. This is quantified by the ratio between the HARP-CDP data and 
the global interpolation shown in Fig. 
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Figure 13: Ratio between the interpolated invariant 7r+ cross sections from HARP-CDP ['41 
and the global interpolation as a function of piab for the angular range 25 < 6'iab < 1 12 degrees 



These deviations are rather consistent with the ones found for protons. Also in this case 
a rapid variation of the cross sections with increasing beam momentum can be excluded by the 
comparison with the pion data from the Serpukhov experiment [iT3]| between 17 and 67 GeV/c 
beam momentum, see Sect. 19. 2 1 below. 

6.2 cos(^^iab) dependence 

As already shown in Sect. 15.21 for protons, the angular distributions at I/a/s = 0.05, in 
between the Fermilab [[TTIl and NA49 |[T2l energies, are presented in Fig. [TH This allows the 
comparison of the two data sets and their connection across the gap in lab angles between 40 
and 70 degrees which represent the upper and lower limit of the respective experiment. 

Further angular distributions at four 1/ \fs values between 0.1 and 0.4 GeV^^ are given 
in Fig. [121 

The angular distributions are characterized by a smooth, close to exponential shape. At 
backward angles, the piab dependence is very steep with four orders of magnitude already be- 
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COS (0|3b) 

Figure 14: Invariant 7r+ cross sections at 1/ -/s = 0.05 as a function of cos(^^iab) combining the 
Fermilab [1 IJ and NA49 [12J data for piab between 0.2 and 1.2 GeV/c. The global interpolation 
is shown as full lines. The measured cross sections in the angular ranges from 70 to 160 degrees 
( BID) and from 10 to 40 degrees ( fl2^ ) are given on the vertical broken lines 

tween piab = 0.2 and 0.8 GeV/c. In forward direction this dependence is much reduced with less 
than one order of magnitude between piab = 0.2 and 1 .2 GeV/c. This is due to the prevailance of 
target fragmentation in this region, see Sect.[10]for a quantitative study of this phenomenology. 

7 The data for negative pions 

This section follows closely the discussion of the 7r+ cross sections in the preceding 
section. 

7.1 1 / dependence 

The invariant tt" cross sections are shown in Fig.[T6]as a function of 1/ y/s for the standard 
grid of ten lab angles between 25 and 1 80 degrees and for constant lab momenta between 0.2 and 

1.2 GeV/c. The interpolated data points in each panel are identified by symbols corresponding 
to the different experiments. 

The solid lines represent the global interpolation by eyeball fits to the data, with several 
features which are worth noticing: 
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COS (0, J COS (0, J 

Figure 15: Invariant 7r+ cross sections as a function of cos(0iab) for four values of 1/ yfs: a) 0. 1, 
b) 0.2, c) 0.3, d) 0.4 GeV~^ and for piab values between 0.2 and 1.2 GeV/c. The standard grid 
of 10 angles, Fig.|71 is indicated by the vertical broken lines 

- All the different data sets form a consistent ensemble without the systematic deviations 
visible in some regions of the proton and 7r+ results. 

- The approach to large beam momenta happens from below for all p\ah- 

- The s-dependence is in general stronger than for 7r+, Fig. [121 If it is again flat up to 
1 / ^/s ~ 0.2 at low piab, it becomes more pronounced both towards higher piab and in the 
approach to the production threshold at large 1 / -y/i indicating a marked increase of the 
7r+/7r^ ratio. 

- This effect has as physics origin the progressive change of the production mechanism 
from pion exchange at low energy to gluon or Pomeron exchange at SPS energy. This 
will be discussed in relation to the charge ratios in Sect. [H 

It is again interesting to compare the energy dependence to the one observed in p+p 
interactions as presented in Fig.fTTl 

Although for both reactions the asymptotic high energy region is approached from below. 
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Figure 16: Invariant cross sections for vr^ in p+C collisions as a function of 1/ ^/s at fixed piab 
and 9iab. The interpolated data points are indicated by symbols corresponding to the respective 
experiments in each panel. The solid lines represent the global data interpolation 



this comparison shows a stronger s-dependence, at the same lab angle, in p+C than in p+p 
collisions. This is due to the component of nuclear cascading which contributes, in the given 
angular range, with equal strength than the target fragmentation to the total yield (see Sect. [TOl 
below). 

7.2 cos(6'iab) dependence 

As for protons and 7r+ in Figs. [10] and [l4l the tt^ cross sections from the Fermilab ifTTTl 
and NA49 lfT2l experiments are compared and combined as a function of cos(6'iab) in Fig.fTSi 

Further angular distributions at four l/y/s values between 0.1 and 0.4 GeV^^ are given 
in Fig. [IE 

Concerning smoothness and continuity these distributions are similar to the 7r+ data, in- 
cluding the large asymmetry between the forward and backward directions. The reduction of 
the cross sections for vr^ with respect to 7r+ with increasing l/y/s is however very apparent. 
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Figure 17: Invariant tt^ cross sections as a function of 1/ y/s for p+p interactions at the two lab 
angles a) 45, and b) 55 degrees for piab values from 0.2 to 1.2 GeV/c. The data are interpolated 
from Blobel ETl and NA49 H2H . The lines are drawn to guide the eye 
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Figure 18: Invariant tt^ cross sections at l/y/s= 0.05 as a function of cos(6'iab) combining the 
Fermilab [1 IJ and NA49 [1121 data for pi^b between 0.2 and 1.2 GeV/c. The global interpolation 
is shown as full lines. The measured cross sections in the angular ranges from 70 to 160 degrees 
( ifTTTl ) and from 10 to 40 degrees ( [11211 ) are given on the vertical broken lines 
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Figure 19: Invariant vr" cross sections as a function of cos(0iab) for four values of 1/ ^/s■. a) 0. 1, 
b) 0.2, c) 0.3, d) 0.4 GeV~^ and for pi^b values between 0.2 and 1.2 GeV/c. The standard grid 
of 10 angles, Fig.|71 is indicated by the vertical broken lines 



This will be quantified in the following Section onn^ /n ratios. 



8 The 7r+/7r ratio 

As already evoked in Sect. 14.5.3] above, the study of /n^ ratios has two main advan- 
tages. Firstly, in this ratio a major fraction of the experimental systematic uncertainties cancels. 
Secondly, the ratio is constrained by very fundamental and model-independent physics argu- 
ments like charge conservation and isospin symmetry. In addition, its s-dependence is governed 
by the hadronic meson exchange process which leads to a power-law behaviour that will be 
shown to be common to a wide range of interactions. In the following argumentation the ratio 
between the global data interpolation for 7r+ and vr" as described in the preceding Sects. [6] and 
|7]will be used: 
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^±(i/yi,piab,M = 7r-T 

) 

As a by-product, the fluctuation of this ratio as a function of angle and interaction energy 
will allow for the estimation of the local precision of the interpolation procedure. 



8.1 The high energy Umit 

It has been established by numerous experimental results that at collision energies in 
the SPS/Fermilab range and above the hadronic interactions are characterized by the absence 
of charge and flavour exchange. It has also been shown that the feed-over of pions from the 
projectile hemisphere into the backward region of xj? is sharply limited to the range of a;^ ^ - 
0.05, see [|28il for a detailed discussion. This range is outside the coverage in Giab and piab 
considered in this publication. 

It is therefore to be expected that the backward production of pions off an isoscalar nu- 
cleus should be charge-symmetric at high energy. This is indeed verified by the results on pion 
production shown in the preceding sections. It is quantified in Fig. [19] which shows the /n^ 
ratio at 1/ y/s = 0.04 or 330 GeV/c beam momentum for all lab angles and lab momenta treated 
in this publication. This number distribution has a mean value of 1.0125 with an rms deviation 
of 3.2%. This rms value may be seen as a first estimate of the local precision of the three- 
dimensional interpolation scheme at this energy which has been established independently for 
both pion charges. 




Figure 20: tt+Zti- ratio at l/^/s = 0.04 GeV^^ for 25 < Oiab < 162 degrees and 0.2 < 
piab < 1.2 GeV/c 



8.2 Energy, momentum and angle dependence of R± 

With decreasing interaction energy or increasing 1 / y/s the 7r+/7r^ ratio develops a strong 
increase at all lab momenta and lab angles. This is shown in Fig. [2T] which gives R± as a function 
of 1 / y/s for four lab momenta. The ratio of the global data interpolation is given in steps of 0.02 
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in 1/ y/s. At each value of 1/ y/s the number of points corresponds to the standard grid of angles 
available at this energy. 
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Figure 21 : R± as a function of 1/ y/s for four values of p\ab- The dots represent the ratio of the 
global data interpolation for 7r+ and tt^ in steps of 0.02 in 1/y/s and for the chosen grid of 10 
angles, a) piab = 0.2, b) pi^b = 0.4, c) piab = 0.6, d) piab = 0.8 GeV/c 



Several features of Fig. [2T] are noteworthy: 

- Considering the wide range of lab angles, R± is at each value of l/y/s confined to a 
narrow band indicating an approximative angle independence. 

- Large R± values in excess of 5 are reached at the upper limit of the available scale in 

- There is a systematic increase of R± with piab. 



8.2.1 Mean tt^/tt ratios and estimation of the local systematic fluctuations of the 
interpolation process 

The features pointed out above may be quantified and at the same time the local system- 
atic fluctuations of the interpolation may be estimated by establishing the mean values {R±) 
averaged over the angular range at each 1 / y/s. These mean values are well defined as shown 
in Fig. [22] which presents the normalized distribution of the point-by-point deviations from the 
mean in percent. 
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Figure 22: Number distribution of the normalized percent deviations AR± from the mean over 
the angular range, a) l/y/s = 0.1, b) l/v/i = 0.2, c) l/i/i = 0.3, d)l/^/s = 0.4. The four piab 
values shown in Fig. [2T|are summed up 



These distributions are of Gaussian shape with an rms which increases with 1 / ^/s as 
indicated in Fig.[ 
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Figure 23: Rms values of the number distributions of the normalized point-by-point deviations 
from the mean charge ratio {R±) as a function of l/\/s. The full line represents a hand in- 
terpolation. Broken line: corresponding error margin of the mean values {R±). Dotted lines: 
corresponding errors for 7r+ and vr^ separately 
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The observed energy dependence of the rms deviations is due to the fact that the invariant 
pion cross sections decrease, after a relatively flat behaviour up to 1/ ^/s ~ 0.15, progressively 
steeper towards the production threshold, see Figs. [121 and [161 This leads inevitably to larger 
variations in the corresponding energy interpolation. 

From the rms values given in Fig. [23lthe error of {R±) may be derived which varies 
between 1% and 5% for the highest and lowest interaction energy, respectively (broken line). 
Also the corresponding error margins for the mean pion yields may be extracted as indicated 
by the dotted lines in Fig. [231 From these plots it appears that the global interpolation induces 
fluctuations which increase from a few percent in the high s region to about 10% in the approach 
to the pion threshold. 

8.2.2 Dependence of R± on Qiab 

The dependence of R± on 6iab is shown in Fig. [24l for four values of 1/ y/s and four 
values of piab- 
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Figure 24: R± as a function of 6iab for four values of l/\/s, a) piab = 0.2, b) piab = 0.4, c) 
Piab = 0.6, and d) pi^h = 0.8 GeV/c. The mean values {R±) are indicated as the horizontal broken 
lines in each panel 
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Evidently no systematic 6iab dependence is visible over the complete angular range 
within the quoted errors. 

8.2.3 Dependence of {R±) on I/a/s andpiab 

In the absence of angular dependence of R± as shown above, the mean values {R±) may 
now be used in order to establish a precise view of the l/y/s dependence for different piab 
values. This dependence is presented in Fig. [251 




1/\ls[GeV''] 

Figure 25: {R±) as a function of l/\/s for five values of piab between 0.2 and 1.0 GeV/c. The 
full lines are hand interpolations through the data points 

Within the errors of {R±) extracted above, a clear piab dependence is evident superposing 
itself to the strong common increase of {R±) with 1/a/s. This increase may be parametrized 
up to l/y/s ~ 0.3 by the functional from 1 + c/s^'^^'^''^ which is, as discussed below, typical 
of meson exchange processes. Indeed the exponent beta varies from 2 to 1.2 for piab increasing 
from 0.2 to 0.8 GeV/c. 

8.3 Interpretation of the observed energy and momentum dependences 

The strong increase of R± with 1 / ^/s merits a detailed study as it is directly connected 
to the basic hadronic production mechanisms in p+A interactions. The fact that the pion yields 
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in the complete backward fragmentation region of an isoscalar nucleus remember the isospin of 
the projectile is clearly incompatible with charge and flavour independent exchange processes. 
Instead a meson exchange mechanism may be invoked which has indeed been used successfully 
in a wide range of work at low projectile momenta, see for instance [18] and references therein. 
Close to the pion production threshold in the nuclear hemisphere, single excitation processes 
via pion exchange of the type 
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vr° (11) 
7r+ (12) 

only allow 7r+ and 7r° production, whereas tt^ production needs double excitation like 
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with in general an additional penalty for tt due to the isospin Clebsch-Gordan coefficients. 
All meson exchange mechanisms are characterized by a strong decrease with projectile energy. 
This energy dependence and its interplay with processes governing the high energy sector is 
studied here for the first time in p+A collisions using the 7r+/7r~ ratio. 

In this context it seems mandatory to first refer to the study of exclusive charge exchange 
reactions in elementary nucleon-nucleon collisions as the complete energy range discussed here 
has been covered there by a number of experiments Il29] - l35l . 



8.3.1 The charge exchange mechanism in elementary nucleon-nucleon collisions 

Charge exchange processes may be cleanly isolated experimentally in nucleon-nucleon 
interactions by studying the following exclusive channels: 

- Charge exchange scattering of the elastic type 

n + p -> p + n (16) 

- Single dissociation with pion production 

p + p n + A++ n+(p + 7r+) (17) 

- Double dissociation with pion production 

P + P (p + TT") + (p + 7r+) (18) 

These channels are characterized by a very steep energy dependence. 
This is to be confronted with non-charge-exchange exclusive channels like: 

- Elastic scattering 

p + p ^ p + p (19) 

- Single dissociation 

p + p ^ p + (p + 7r+ + TT") (20) 
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- Double dissociation 



p + p (p + 7r+ + TT ) + (p + 7r+ + tt" 



(21) 



which show a constant or logarithmically increasing s-dependence. 

Charge exchange scattering has been measured by five experiments in the range of neu- 
tron beam momenta from 3 to 300 GeV/c. ||29] - [33l . This is exactly covering the energy range 
discussed in this paper. The single and double dissociation has been studied at the CERN ISR 
by two experiments [|34l[35l extending the energy scale to s = 3700 GeV^. The two ISR ex- 
periments may be directly compared to the charge exchange measurements after appropriate 
re-normalization of the cross sections in the overlap region at the lowest ISR energy. 

The resulting s-dependence at a momentum transfer t = 0.032 GeV^ is presented in 

Fig.r 
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Figure 26: Invariant cross sections of charge exchange and single and double dissociation in 
nucleon-nucleon interactions as a function of s at a momentum transfer t = 0.032 GeV^. The 
full line represents an interpolation of the data points. The insert gives the local slope (3 in the 
parametrization / ~ s^^ as a function of s 
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Several features of Fig. [26] are of interest: 

- There is a decrease of about 4 orders of magnitude in cross section between the lowest 
and highest s value. This decrease is to be compared to the constant or logarithmically in- 
creasing elastic and nucleon diffraction cross sections. The charge exchange contribution 
is therefore negligible compared to the inclusive baryon yields already at SPS energy. 

- There is a steady decrease of the local slope df / ds with energy, from about 3.6 at 3 GeV/c 
to about 1 . 1 above 80 GeV/c beam momentum. 

- A characteristic change of slope manifests itself at around 30 GeV/c beam momentum. 
These features have been interpreted in the 1970's when the relevant experiments were 

performed, in the framework of Regge theory which predicts an s-dependence of the form 

/ ^ = s-^, (22) 

where a is the intercept of the leading trajectory. This should in the case of one-meson exchange 
at low energy be given by the pion trajectory with zero intercept. The actual beta values above 
3 at low s seem to contradict however this expectation. Here threshold effects may play a role 
which are not included in the parametrization ((22l) . 

With increasing energy the slopes move through the region of pion exchange with /3 ~ 2 
down to values of about 1.1 at high energy which could be connected to p and a2 exchange 
with correspondingly higher intercepts a in the region of 0.5. At ISR energy the ratio of p/n 
contributions has indeed been estimated to be about 2 [35 1. Anyway the simple parametrization 
given by ([22l) should not be expected to hold over the full energy scale. What is interesting 
here is rather the strong decline of the charge exchange cross sections with energy and the 
experimentally rather precisely determined slope variation. 

8.3.2 A remark concerning baryon resonance production in hadronic interactions 

The single (flTI) and double (fTSi) dissociation processes defined above are determined by 
the formation of A resonances in the final states. They therefore constitute a source of direct 
A production in nucleon-nucleon interactions. These channel cross sections decrease rapidly to 
the fxbam level at SPS energies. In contrast, the non-charge exchange channels like ([201) and 
([2T[) have no s-dependence and stay on the mb level of cross sections. Their final states have 
been shown to be governed by N* resonances [|36l which may be excited by Pomeron exchange. 
Moreover, the p-i-7r+ combination of the p-i-7r+-i-7r" final states has been shown to be dominated 
by A++ [[371 . This is an indirect source of A resonances as a decay product of N* states which 
have large decay branching fractions into A-i-tt and A+p. It is therefore questionable if, at SPS 
energies and above, any direct A production is persisting. This is an interesting question for 
the majority of microscopic models which produce final states by string fragmentation. In the 
baryonic sector, diquark fragmentation is generally invoked with a prevailing direct production 
of A resonances which by isospin counting will dominate over N*. Indeed in practically all 
such models there is no or only negligible N* production. As shown below, the decrease of 
charge exchange processes can be traced well into the non-diffractive, inelastic region of particle 
production. The multi-step, cascading decay of primordial N* resonances into A resonances and 
final state baryons should therefore be seriously considered, in particular also concerning the 
consequences for the evolution of final state energy densities with time. 

8.3.3 The charge exchange mechanism inp+C interactions as a function of interaction energy 

The very characteristic decrease of {R±) with increasing s derived from the global data 
interpolation. Fig. [25l offers a tempting possibility of comparison to the phenomenology dis- 
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cussed above for the elementary nucleon-nucleon sector. Indeed, two components should con- 
tribute to the observed /n^ ratios: at high energy this ratio should approach unity due to the 
absence of charge and flavour exchange in this region. At low energy on the contrary it should 
be governed by meson exchange with its strong s-dependence. These two components may be 
tentatively separated by using instead of {R±) the quantity 

{RT) = {R±) - 1 (23) 

in order to extract the meson exchange contribution. This quantity is plotted in Fig. [27] as a 
function of s for four pi^b values from 0.2 to 0.8 GeV/c. 
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Figure 27: (R^^) as a function of s for different values of piab, a) pi^h = 0.2, b) piab = 0.4, c) 
Piab = 0.6, and d) piab = 0.8 GeV/c. The elastic limit is indicated by the arrows 

A very characteristic pattern emerges which resembles the s-dependence for the charge 
exchange in elementary interactions described above, see Fig. [26l In general (i?™*^) follows a 
power law dependence on s 

{R"^') ~ cs"^™"^ (24) 

with local slopes (3^^ which are in turn a function of s. Three different regions with distinct 
local slopes can be identified in Fig. [271 

- A first region with large slopes is located at s below about 6 GeV^. This region is strongly 
influenced by threshold effects as the threshold for inelastic production is placed at the 
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elastic limit s = Ani^ = 3.5GeV^ indicated in Fig.[27l In the approach to pion thresh- 
old the TT^/n^ ratio has to diverge as tt" is progressively suppressed, see above. With 
increasing pi^b this suppression will of course be more pronounced. 

- An intermediate region between about 8 and 40 GeV^ with an s dependence decreasing 
with increasing pi^^. 

- A third region with flattening s-dependence above about 40 GeV^. 

At the lowest piab value of 0.2 GeV/c corresponding to the lowest momentum transfer, the 
similarity to the charge exchange process in nucleon-nucleon interactions. Fig. [261 is absolutely 
striking. This concerns both the detailed shape and the overall suppression factors. With in- 
creasing piab, the s dependence is modified in a systematic way by a general reduction of slopes, 
with the exception of the threshold enhancement. This is quantified in Fig. [28] which shows the 
local slopes as a function of s for pi^b values between 0.2 and 1 GeV/c. 
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Figure 28: Slopes of the s-dependence of (i?™'^) as a function of s: a) p^b = 0.2, b) p^b = 0.4, 
c) piab = 0.6, d) piab = 0.8, and e) p^b = 1.0 GeV/c. The shaded regions mark the error margins 

With the exception of the threshold region, the slopes are confined to the region between 
2 and 1 typical of meson exchange processes. The dependence on pi^b is given in Fig. [29] where 
the slopes in the three regions of s specified above are presented. 

This Figure shows clearly the different nature of the low s enhancement where the slopes 
increase strongly with pi^b- The two other regions, full and dotted lines, are compatible with 
a Regge parametrization with trajectory intercepts which increase with piab. This is insofar in- 
teresting as the region of measurements regarded here covers the complete backward angular 
range and the corresponding interactions are by no means confined to diffractive or low momen- 
tum transfer collisions. It is shown in Sect. [10] of this paper that in the backward hemisphere 
the pion yields from nuclear cascading and target fragmentation are comparable. If the nuclear 
component is characterized by low momentum transfer reactions [[3 the target fragmentation 
is manifestly inelastic and non-diffractive. It governs the total yield at all angles below about 
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Figure 29: Local slopes in the regions s < 6 GeV^ (broken line), 6 < s < 40 GeV^ (full line) 
and s > 40 GeV^ (dotted line) as a function of piab between 0.2 and 1 GeV/c 



70 degrees. 

In conclusion of this study of 7r+/7r" ratios in p+C interactions the following points 
should be stressed: 

- The global data interpolation leads to a precise and consistent description of the behaviour 
of the vt+Ztt" ratios in the full backward hemisphere, thus offering an additional tool for 
the discrimination of experimental deviations. 

- The inspection of the detailed s-dependence of the ratios opens a new window on the 
underlying exchange processes. 

- In particular the comparison to the elementary nucleon-nucleon collisions establishes a 
close relation between apparently disjoint sectors of the different hadronic interactions. 

9 Data sets not used in the global interpolation 

As mentioned in Sect. 14. 51 four of the 19 investigated data sets are incompatible with the 
attempt at generating an overall consistent description of the experimental situation. These data 
will be shortly discussed below. 

9.1 The proton data of ref. Q 

These data have been obtained at the Bevalac using beam momenta of 1.75, 2.89 and 
5.89 GeV/c, spanning a lab momentum range from 0.3 to 0.9 GeV/c at a lab angle of 180 de- 
grees. The resulting cross sections trace the shape of the 1/ \fs dependence rather precisely but 
are consistently about a factor of two below the global interpolation as shown in Fig. [30l Here 
the full lines correspond to the global interpolation and the broken lines give the interpolation 
divided by a factor of two. 

As the angular bin from 160 to 180 degrees is mostly covered by data around 160-162 de- 
grees, a steep angular dependence in this region cannot a priori be excluded. The smooth and 
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Figure 30: The proton cross sections from fj] in comparison to the global data interpolation (full 
lines) in the lab angle bin from 160 to 180 degrees. The broken lines correspond to a reduction 
of the interpolation by a factor of two. The inserted histogram gives the number distribution of 
the ratio between data and interpolation 



gentle angular dependence of the interpolated data shown in Fig. [3T] for the angular range from 
82 to 180 degrees and for the three l/y/s values of ref. [|7||, together with the constraint of the 
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Figure 3 1 : Invariant proton cross sections as a function of cos(9iab) for three values of 1/ y/s; a) 
0.427, b) 0.366 and c) 0.281 GeV^^. The full lines give the global data interpolation, the open 
circles the data from tVJ 
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approach to 180 degrees with tangent zero, excludes however a drop of the cross sections by a 
factor of two between 160 and 180 degrees. 

9.2 The proton and pion data of refs. ll3l[T3]l 

A sizeable set of data on proton |l3l and pion lfT3ll production has been obtained at the 
Serpukhov accelerator spanning the range of beam momenta between 17 and 57 GeV/c. This 
fills the gap between the PS and SPS energies where no other data are available. The data cover 
the piab range from 0.25 to 1.2 GeV/c at 6iab = 159 degrees. They are presented in Fig. [32] in 
comparison to the global data interpolation. 




1 /\ls [GeV^] 1 /\|s [GeV^] 1 /\|s [GeV^] 

Figure 32: The data of [I3l[l3l as a function of at Giab = 159 degrees (open circles) in 

comparison to the global data interpolation at 160 degrees (full lines) 

Several features of this comparison are noteworthy: 

- The shape of the dependences complies precisely with the global interpolation. 
This is compatible with the absence of rapid variations of the cross sections with energy 
in the region between PS and SPS. 

- There is a pronounced suppression of these data with respect to the interpolation with 
increasing pi^b reaching factors of three at the upper ranges for protons and pions. 

- The 7r+ and tt^ data show an identical behaviour. 

- The proton data are tracing the interpolation up to piab = 0.4 GeV/c whereas the pion data 
are already suppressed in this piab range. 

- The suppression factors are generally bigger for the pions at equal pi^b- 

- These features might be compatible with a momentum scale error. 

In addition to the reproduction of the shape of the l/\/s dependence, also the n^/n^ 
ratio complies exactly with the one extracted from the global interpolation. Fig. [331 up to 
Piab = 0.7 GeV/c. Above this value there is a sharp drop of R± reaching unphysical values 
at the upper limit of piab. This drop of about 20% has however to be compared to a drop of 
300% of the invariant cross sections at this limit. 
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Figure 33: Pion charge ratio R± from |fT3l as a function of piab averaged over beam momentum 
between 17 and 57 GeV/c. The full line gives the result of the global interpolation averaged 
over the same beam momentum scale, Sect. [8] 



9.3 The pion data of ref. lH 

These results cover a range from 3 to 12 GeV/c beam momentum at ©lab between 25 and 
117 degrees and 0.125 < pi^b < 0.75 GeV/c. They are thus directly comparable to the ones 
from [31 which are part of the global data interpolation. Their differences to this interpolation 
are presented in Fig. [34]for all beam momenta and the standard grid of p^h and 6iab values. 
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Figure 34: Histograms of the percent differences for all angles and beam momenta between Bl 
and [|15(l . Panel a) tt^, panel b) tt"*" 
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If the mean values of the differences are close to zero, their number distributions show 
wide spreads especially for 7r+. This is exemplified in Fig. [35] where a typical comparison to the 
global interpolation (full lines) is given as a function of 1/ y/s at ©lab = 67 degrees for four piab 
values. 
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Figure 35: Invariant pion cross sections from ifTSl for 6iab = 67 degrees and four pi^b values as 
a function of l/y/s, (open circles) in comparison to the global data interpolation (full lines). 
Panel a) 7r+, panel b) tt^ 
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A comparison of tx^ jix ratios is given in Fig. 

©lab atpbeam = 5 GcV/C. 



as a function of piab for four values of 
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Figure 36: Charge ratio R± from [1151 as a function of piab for four lab angles and a beam 
momentum of 5 GeV/c. The full lines give the results of the data interpolation from [15 J, the 
full circles correspond to the ratios of the measured cross sections. The broken lines give the 
result of the global interpolation 



This figure demonstrates the importance of using, in addition to the invariant cross section 
proper, the particle ratios which are strongly constrained by physical arguments, see Sect.dl 

9.4 The pion data of ref. BH 

These data have been obtained at a beam momentum of 3 1 GeV/c in a 6iab range from 
0.6 to 22.3 degrees and piab from 0.2 to 18 GeV/c. If a large part of the given angular and 
momentum coverage falls outside the backward region regarded here, the low momentum range 
up to piab ~ 0.5 GeV/c for all angles and the range 0.6 < piab < 1 GeV/c for angles above about 
9 degrees corresponds to negative xp and can therefore be considered here. 

The complete 1 / ^Js dependence established in the preceding sections has a lower angular 
limit at 25 degrees corresponding to the lowest value of the standard grid of angles. This angle 
is close to the highest angle of [fT4l at 22.3 degrees allowing for a safe interpolation. This is 
shown in Fig. [37] where the global interpolation is compared to the data of lfT4l at their highest 
angles between 12 and 22 degrees for two piab values both for 7r+ and for vr". 

As the global interpolation is limited to ©lab > 25 degrees, another way of comparison is 
offered by the combined NA49 and Fermilab results at 158/400 GeV/c where the former data 
cover the complete angular range of [[T4]| . The ratio of the available high energy data to the 
results at 3 1 GeV/c beam momentum is shown, as a function of 6iab, in Fig. [381 
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Figure 37: Invariant pion cross sections as a function of cos(6iab) at l/\/s = 0.13 GeV^^ in 
the 6iab region from 12 to 35 degrees for piab = 0.4 and 0.6 GeV/c. Full lines: global data 
interpolation. Full dots and broken lines: data from lfT4ll 
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Figure 38: Ratio of invariant pion cross sections at 158/400 and 31 GeV/c as a function of 
©lab for Piab = 0-4 and 0.6 GeV/c. The full circles correspond to the global data interpolation, 
the open circles to the direct ratio ^M- The broken lines give the mean ratio over the 
complete range of the global data interpolation 
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Apparently this cross section ratio is within errors angle independent over the full range 
of the global data survey, with well defined averages below 1.05 for 7r+ and 1.1 for tt^. In 
contrast, the cross section ratio between NA49 and ref. [fT4ll shows values in the region of 1 .4 
increasing with decreasing ©lab. 

In conclusion to this Section it may be stated that the global data interpolation between 
15 different experiments attempted in this paper proves to be a useful tool for the detection 
of deviating data sets. Further details concerning the above comparisons can be found in an 
internal report on ref. [2611. 

10 The separation of target fragmentation and intra-nuclear component for pion 
production at SPS energy 

Hadronic production in the backward direction of p+A collisions has two components: 
the fragmentation of the target nucleons which have been hit by the projectile proton, and the 
propagation of momentum transfer into the nucleus by secondary nucleon-nucleon interaction 
which follow, on a longer time scale, the initial excitation process. Both processes are governed 
by the mean number of collisions (z/) suffered by the projectile on his trajectory through the 
nucleus. 

As only the sum of these two separate mechanisms is experimentally accessible, a min- 
imum assumption about the fragmentation of the target nucleons is needed in order to allow 
the separation of the components in an otherwise model-independent fashion. This minimal as- 
sumption consists in assuming that the fragmentation process of the hit nucleons is equal to the 
basic nucleon-nucleon interaction, taking full account of course of isospin symmetry. In addi- 
tion and only valid for the relatively small value of (y) in the Carbon nucleus, it will be assumed 
that successive collisions result in hadronization at full interaction energy of the corresponding 
elementary interactions. 

As far as the value of (z/) is concerned, this has been determined for pion production in 
some detail in [28| using the forward and the backward region dX xp > -0.1 where no intra- 
nuclear cascading is present, see below. This determination used three independent approaches: 

- A Monte-Carlo calculation using the measured nuclear density distributions. 

- The relation between the inelastic cross sections of p-i-p and p+C interactions. 

- The approach to x^? = -0.1 of the ratio of pion densities in p+C and p-i-p collisions. 

The two former methods have to make the assumption that the inelastic interaction cross 
sections are independent of the number of subsequent collisions u. 

In [|2l| a similar approach is used concerning the production of protons and anti-protons, 
again in the regions where there is no contribution from nuclear cascading as well as in the full 
backward hemisphere. 

All methods mentioned above result in a consistent estimate of (z/) = 1.6 in p+C colli- 
sions, with a relative systematic uncertainty of the order of a few percent. 

In the following argumentation a prediction of the mean pion density of target fragmen- 
tation in the backward hemisphere at ^ = 17.2 GeV will be used which is relying on the 
published pion data from NA49 [21] and the estimated mean number of collisions, (z/). The 
invariant pion cross sections are divided by the inelastic cross section to yield the quantity 

{U{xf,Pt)) = 0.5{f;^\xF,pT) + fl{xF,PT)) (25) 
per inelastic event which establishes isospin symmetry, and 

r''{xp,pT) = l.6{U,{xF,PT)) (26) 



42 



This prediction is transformed into the appropriate coordinates pi^b and 6iab and divided 
by the measured invariant p+C cross sections /pc(piab, ©lab) per inelastic event yielding the ratio 



-K [Pub, flab J — 7^; T 

JpCyPlab, c'labj 



(27) 



This ratio is shown in Fig. [39|as a function of piab for the lab angles 10, 20, 30, 40 and 
45 degrees. 
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Figure 39: R^'^^'^ipiab, ©lab) as a function of piab for the five angles 10, 20, 30, 40 and 45 degrees 



It is evident that the ratio is close to one for the three lowest angles at all pi^b and for the 
region below 0.8 GeV/c for 40 and 45 degrees. This is quantified in Fig. l40l which gives the 
distribution of the ratio for the mentioned piab ranges. 

The results show that indeed the measured pion cross sections correspond for lab angles 
up to 45 degrees precisely to the prediction from elementary collisions. This indicates that 
there is no contribution from intra-nuclear cascading in this region, in accordance with the 
results of [|28l . A drop of the ratio becomes however visible in the higher piab range at 40 and 
45 degrees. This marks the onset of a nuclear component which becomes clearly visible in the 
ratios at larger angles shown in Fig.l4Tl 

It is interesting to note that the target fragmentation governs the pion density up to the 
highest lab angles at low piab, with RP^^'^ values of more than 50%. The ratio decreases how- 
ever steadily with increasing piab and reaches zero at distinct momentum values indicating the 
approach to = -1 in the plots of Fig. [TJ This kinematic effect is more clearly brought out 
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Figure 40: Distribution of RP^^'^{piab,Qiab) for the lab angles 10, 20, 30 degrees and 40 and 
45 degrees at piab < 0.8 GeV/c. Full line: Gauss fit to the distribution with a mean value at 1.01 
and a relative rms of 4.3% 
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Figure 41: -R'"''^'^(piab, ©lab) as a function of pi^b for the angles of 45, 55, 70, 90, 112, 137 and 
160 degrees. The full lines are local interpolations 



in Fig. |42] showing that the fraction of target fragmentation is essentially a function of and 
rather independent on lab angle. 
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Figure 42: a) measured correlation between piab and Giab for constant values of R'^'^^'^ipuh, ©lab) 
between 0.1 and 0.9, b) correlation between pi^b and 6iab for fixed values of xp between -0.1 
and -1.0 and c) BP'^^'^ipiab, ©lab) as a function of xp 

The correlation between pi^b and 9iab for fixed values of RP^'^'^ shown in panel a traces 
rather exactly the kinematic correlation between the same variables for fixed values of xp, 
panel b. This allows to establish a direct dependence of EF'^^'^ on xp which is to first order 
angle-independent, panel c. 

The invariant densities P'^^'^{p\ab, ©lab) per inelastic event as predicted from the fragmen- 
tation of the participant target nucleons is presented in Fig. |43l 




COS(0|ai,) 

Figure 43: Predicted invariant density /'"^'''^(piab, ©lab) per inelastic event as a function of 
cos (6iab) for fixed values of piab between 0.1 and 1.2 GeV/c. The full lines represent data inter- 
polations 
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This density may be subtracted from the pion density /(piab, 0iab)/cr'"'^' measured in p+C 
interactions which is within errors equal for 7r+ and vr^, see Figs. [14]and[T8l The resulting 
invariant density 

n\pu,, e,ab) = ^^^^fr^ - r^'iPU^, ©lab) (28) 

is shown in Fig. |44l 




This subtraction procedure becomes of course uncertain in the small angle region where 
the nuclear component is on the few percent level and below with respect to the target fragmen- 
tation, see Figs. [39] and l4n 

The invariant angular distributions shown in Figs. I43]andl44] may be converted into num- 
ber distributions following: 

1 T^; = 27r— /'^ (plab, Blab) (29) 

apiabttfe'iab -C/lab 



and 



(i^n""'^'(piab, ©lab) r,^P\ab nnuclf^ o N 

1 = 2vr — / (Plab, ©lab) (30) 

apiabttfe'iab -C/lab 

Integrating these distributions over piab, the number distributions 
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and 



(icos(9ia 



(31) 



(i cos (Slab) 

are obtained which are shown in Fig. |45] together with the ratio 



i?™'^\cos(eiab)) 



(icos(6iab)/ (icos(9iab) 
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Figure 45: a) dn"^^^ / d cos{<diab) as a function of cos(6iab) (full line), (inP'''^''/c/cos(0iab) as a 
function of cos(6iab) (broken line), b) the ratio i?""'^' as a function of cos(9iab) 

Evidently the nuclear component of pion production stays comparable to the target frag- 
mentation in the full backward hemisphere of Oiab- It decreases rapidly for ©lab below about 
60 degrees and vanishes below ©lab 25 degrees. 

Integration of (in™'^'/(icos(0iab) over cos(6iab) results in the total single pion yield from 
nuclear cascading 



= 0.105 



per inelastic event. The predicted integrated yield from target fragmentation is 



pred 



with 



and 



4 



n^l = 3.018 



2.151 



(34) 



(35) 



(36) 



= 2.360 



(37) 
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from p+p interactions as measured by NA49, [ED. This means that for p+C interactions the 
nuclear component of pion production amounts to 4.9% of the pions originating from the frag- 
mentation of the hit target nucleons. Applying isospin symmetry on the isoscalar C nucleus 
with 



?7,^+ = n^- = rij^o (38) 

the total pion yields are 6.45 from target fragmentation and 0.315 from nuclear cascading. 

Making use of the kinematic relation between the coordinate pairs piab, ©lab and xp, Pt, 
see Fig.[l]:, the double differential yields for the nuclear component as functions of Xi? and px 

-27rp^..^n\xF,pT) (39) 



dxpdpT E 

may be obtained where Pmax, ®, and E are cms quantities. The resulting pion density distribu- 
tions are shown in Fig.|46]as a function of xp for pt values from 0.05 to 0.7 GeV/c. 




Figure 46: Double differential pion density cPn"-^'^^ / dxpdpT as a function of a;^ for a) 0.05 < 
Pt < 0.3 GeV/c and b) 0.4 < pr < 0.7 GeV/c 



A peak at low px and xp = -0.15 is apparent which corresponds to the location of pi- 
ons with small lab momentum, see Fig.[T] With increasing pt the maximum density decreases 
and shifts in to lower values which is again in accordance with the kinematic correlation 
visible in Fig. \T\ Integration over px results in the single differential density dn°^'^^ /dxp^xp) 
shown in Fig. |47] together with the predicted density distribution dnP^^'^ /dxp{xp) from target 
fragmentation and with the ratio of the two densities. 

The Pt integrated pion density dn^^'^^ /dxp{xp) shows a peak at xp ~ -0.2 and van- 
ishes at Xp -0.08. As shown by the density ratio with the predicted target fragmentation 
dnJ"^^'^/dxp{xp) in Fig.l47b. the nuclear component reaches 10% of the target fragmentation at 
Xp = -0.15 and exceeds this contribution for x^? < -0.55. 

The nuclear pion component extracted above is used in [|2l in conjunction with the com- 
plementary nuclear proton component to obtain the percentage of cascading protons which are 
accompanied by pion emission. 
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Figure 47: a) Pion density drf"^^/dxp as a function of xp (full line). The predicted den- 
sity distribution from target fragmentation dnP^^^/dxp is shown as the broken line; b) Ratio 

R'^^^^xf) = {dn"'"'^/dxF)/{d'nP"'^/dxF) as a function of 



11 Conclusion 

This paper presents a survey of available data concerning backward proton and pion pro- 
duction in minimum bias p-i-C interactions, including new and extensive data sets obtained at 
the CERN PS and SPS. The backward direction being defined as the complete phase space at 
negative Feynman xf, the data cover, for projectile momenta from 1 to 400 GeV/c, the ranges 
from 0.2 to 1.2 GeV/c in lab momentum pi^b and from 10 to 180 degrees in lab angle Oiab- 
The paper attempts an interconnection of the different data sets by a detailed three-dimensional 
interpolation scheme in the variables l/i/i, piab, and cos(6iab). This attempt allows a precise 
control of the internal data consistency as well as the study of the evolution of the invariant 
inclusive cross sections in all three variables. 

A literature search has provided a set of 19 different experiments with a total of more 
than 3500 data points. These measurements were obtained over 40 years of experimentation 
by collaborations employing widely different experimental techniques. In this respect it may 
be stated as a first positive result that the majority of the data may be combined into a sur- 
prisingly self-consistent ensemble. This global interpolation scheme results in a considerable 
discriminative power against the systematic deviation of particular data sets. Only 4 of the 19 
quoted experiments show in fact deviations which clearly mark them as systematically diverg- 
ing. These experiments are inspected in detail one by one in an attempt to clearly bring out the 
discrepancies. In some of the cases, possible experimental error sources are pointed out. 

The underlying physics provides for additional constraints concerning basic quantities 
like charge conservation and isospin symmetry as well as the necessity of smoothness and 
continuity of the observed cross sections. Whenever possible, contact to the complementary 
elementary nucleon-nucleon interactions is established. This concerns in particular the evoca- 
tion of mesonic exchange processes for the description of /n^ ratios and the prediction of 
the target fragmentation from elementary interactions and its separation from the component of 
nuclear cascading. 

As far as the dependences of the invariant cross sections on the three basic variables piab, 
©lab and is concerned, a well constrained phenomenology emerges. The piab dependences 
are exponential or close to exponential over a major part of the phase space with some excep- 
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tions mostly towards low interaction energies. This fact results in an important constraint for 
the data interpolation. The cos(6iab) dependences are not far from exponential and smooth and 
continuous through all lab angles. In particular there is no indication of an instability around 
90 degrees for the proton yields. The 1 / ^ dependences converge, after strong variations close 
to production threshold, smoothly to asymptotic behaviour in the SPS energy range. This re- 
gion is approached from above by the protons and from below for the pions. This convergence 
is confirmed by the 7r+/7r^ ratios which show, being governed by meson exchange at low -y/i 
with large values marked by the projectile isospin, a smooth decline with energy towards unity 
as expected from the underlying elementary exchange processes. 
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